Raman scattering in a Heisenberg 5 = 1/2 antiferromagnet on the triangular lattice 
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We investigate two-magnon Raman scattering from the 5 = 1/2 Heisenberg antiferromagnet on the triangu- 
lar lattice, considering both the effect of renormalization of the one-magnon spectrum by 1/5 corrections and 
final-state magnon-magnon interactions. The bare Raman intensity displays two peaks related to one-magnon 
van-Hove singularities. We find that 1/5 self-energy corrections to the one-magnon spectrum strongly mod- 
ify this intensity profile. The central Raman-peak is significantly enhanced due to plateaus in the magnon 
dispersion, the high frequency peak is suppressed due to magnon damping, and the overall spectral support 
narrows considerably. Additionally we investigate final-state interactions by solving the Bethe-Salpeter equa- 
tion to 0(1/ S). In contrast to collinear antiferromagnets, the non-collinear nature of the magnetic ground state 
leads to an irreducible magnon scattering which is retarded and non-separable already to lowest order. We show 
that final-state interactions lead to a rather broad Raman-continuum centered around approximately twice the 
'roton' -energy. We also discuss the dependence on the scattering geometry. 
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I. INTRODUCTION 

Raman scattering is an effective tool to study the excitation 
spectrum of magnetic systems since the intensity of the in- 
elastically scattered light is directly related to the density of 
singlet states at zero momentum. In local-moment magnets 
with well defined magnon excitations this quantity is linked 
to the two-magnon density of states. Therefore, magnetic Ra- 
man scattering plays an important role in understanding the 
dynamics and interactions of magnons in conventional spin 



system: 
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This is particularly true for the spin-1/2 square- 



lattice Heisenberg antiferromagnet (HAF) of the high-T c su- 
perconductor parent compounds, where experimental^^ 7 , and 
theoretica l 8,9,10 ' 11 studies of the magnetic correlations by Ra- 
man scattering may provide important insight into the en- 
ergy scales relevant to the pairing mechanism (for reviews see 
Refs^). 

Raman scattering from HAFs can be understood in terms of 
the Loudon-Fleury (LF) processes^, in which two magnons 
are simultaneously created by light absorption and emission. 
In the limit of large on-site Coulomb correlations U, the 
Hamiltonian describing these processes can be obtained as a 
leading term of the expansion in t/(U — u>), where t is the 
nearest-neighbor (NN) hopping, and u> is of the order of pho- 
ton frequencies^. 

The Raman intensity of HAFs on hypercubic lattices with 
unfrustrated NN exchange and collinear type of antiferro- 
magnetic (AFM) order allows for a straightforward semi- 
quantitative interpretation in terms of the LF processes. In 
fact, in real space, exchanging two NN spins of 5 — 1/2 leads 
to an excitation with energy ~ (z — 1) J, where z is the co- 
ordination number and J is the AFM exchange energy. The 
reduction of 0/ J by — 1 is a consequence of the exchange link 
between the NN sites and can be interpreted in terms of a two- 
magnon interactions in the final state. In momentum space, 
the linear spin-wave theory yields non-dispersive magnons 
along the the magnetic Brillouin zone (BZ) boundary, leading 
to a square-root divergence of the bare two-magnon density of 
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states at Q = zJ. Inclusion of the final state magnon-magnon 
interactions broadens the singularity and shifts it down to 

~ 2.9 J 8 ' 9 ' 10 ' 11 in two dimensions, which is consistent, both 
with the real-space result fl = J(4 — 1) = 3 J, and with the 
experimentally observed Raman profile. 

In contrast to conventional collinear HAFs, very little is 
known theoretically about Raman scattering from frustrated 
HAFs. This is intriguing, since the singlet spectrum is be- 
lieved to be an essential fingerprint of such magnets. The spin 
5=1/2 HAF on the triangular lattice (THAF) with NN ex- 
change interactions is a prominent example of strongly frus- 
trated spin systems. It has a ground state with non-collinear 
120° degree ordering of the spins. Due to this non-collinearity 
of the classical ground state, nontrivial corrections to the spin- 
wave spectrum appear already to first order in 1/S*. It has been 
shown in Refs i 16 ' 17,18 that 1/S corrections strongly modify the 
form of the magnon dispersion of the triangular HAF. The re- 
sulting dispersion turns out to be almost flat in a wide range 
of momenta in which it possesses shallow local minima, "ro- 
tons", at the midpoint of the faces of the hexagonal BZ. This 
differs strongly from the classical spin-wave spectrum, which 
lacks such minima and flat zones. Similar results have been 
obtained in series expansion studies 19 . 

Motivated by these recent findings, in this paper, we an- 
alyze the Raman scattering from the THAF by 1/5 expan- 
sion. This is complementary to the recent analysis of Ra- 
man scattering on finite, 16 sites THAFs by means of exact 
diagonalization^ .. First, our results show that the Raman in- 
tensity is very sensitive to both 1 / S corrections of the magnon 
spectrum and the magnon-magnon interactions in the final 
state. Moreover, we find that the Loudon-Fleury process on 
the THAF leads to a Raman profile, which is independent at 
0(1/5) of the scattering geometry. 

The manuscript is organized as follows. In section!]]] we re- 
view results from existing calculation s 16 ' 17 of the one magnon 
excitations in the THAF to first order in 1/5 needed for our 
study of the Raman spectra. In section [Til] we consider the 
LFP to leading order in 1/5. In section HVl we calculate the 
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FIG. 1: a) Classical 120° degree non-collinear spin order on the tri- 
angular lattice. Basic vectors of triangular lattice: S\ = (i -^), 
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and (?3 = (1,0). b) Definition of scattering angles 



for LF vertex. 



Raman spectrum on various levels of approximation in 1/5, 
i.e. bare, one-magnon renormalized, and including final state 
interactions, and show that Raman profile is very sensitive to 
the magnon-magnon interactions. We discuss our results in 
section [V] 



II. MODEL 



The Hamiltonian of the THAF reads 



H 



(ij) 



(1) 



where S, are spin— 1/2 operators, i refers to sites on the trian- 
gular lattice, () denotes NN summation, and J is the exchange 
interaction. The classical ground state of the THAf£i is a 
non-collinear 120° degree ordering of spins which is shown in 
FigOJ)- To avoid the complexity of a three-sublattice notation 
it is convenient to work within a locally rotated frame of refer- 
ence in which the magnetic order is ferromagnetic. To achieve 
this we assume a gauge in which the (x, z)-coordinates la- 
bel the lattice plane and a uniform twist with a pitch vector 
Q = (47r/3, 0) is applied to the y-axis. The laboratory frame- 
of-reference spin Sj is related to the spin S, in the rotated 
frame through 



Si = 







sinfe) - cos(q, 

0-1 
cos(gi) sin(gi) 



(2) 



where = 27r(2/; + mj)/3 and (ij, m;) are integers labeling 
the points on the triangular lattice, which is depicted in Fig. 
QJ). In contrast to Si, the spin Si is amenable to a representa- 
tion in terms of a single Holstein-Primakoff boson field on all 
sites 



S! = S-a+a^ 

5+ = (25 - afa, f' 2 ^ 

S- = a+(2S-a+a t ) 1 ' 2 



(3) 



Because we intend to study magnon interactions to first or- 
der in 1 /5, we need to expand the Hamiltonian in Eqn. up 



to quartic order in the boson fields. We have 

H - E = 3JS{H 2 + Hz + H 4 ) , (4) 
where Eq = 3JS 2 /2 is the classical ground state energy and 



h 2 = Y, A * a W + ^( a k a -k + «k a -k) 



H 4 



(5) 



(6) 



ki ,k 2 ,k 3 



(l4ki + ^k 2 )^k 3 ,ki+k2 



1 

165, 



E *k»+k4,ki+k a 4 I a k B a ks a lu X {T ' 
ki,k2,k 3 ,k 4 

(4(^ kl _k 3 + fk 2 -k 3 ) + ^ki + ^k 2 + ^k 3 +^k 4 ))- 

2<5k 1+ k 2 +k 3 ,k 4 (a kl a k 2 a L a k4 + a k 4 a k 3 a k 2 a k 1 ) 

(z^j + »k 2 +1^3) , 



where momentum k is defined in the first magnetic BZ. We 
use the following notations: 



A k = l + ^ k /2, B k = -3i*/2, 
and the momentum dependent functions are 



(cos k x 



2 cos — cos 

2 2 

ky v^3 

- sin — (cos fez — cos — - 



T 



(8) 

(9) 
(10) 



The expressions for and i?4 have been obtain first in 
Ref»i£. The essential difference between Eqns. - (0 and a 
corresponding expansion around a Neel state on a hypercubic 
lattice is the occurrence of the term H3, which is present due 
to the non-collinearity of the classical ground state configura- 
tion of the THAF. In the remainder of this paper we set the 
scale of energy to 3 J/2 = 1, i.e. for 5=1/2 the prefactor in 
Eqn. (HJi is unity. 

To proceed, we diagonalize the quadratic part of the Hamil- 
tonian H2 by a Bogoliubov transformation to a set of magnon 
quasiparticles 



"kC k + w k c^ k 
"kC k + ^ k c_ k , 



(ID 



where c k are bosons, and the coherence coefficients 



"k 



A k + £ k 



(12) 



2E V 



satisfy u k — v 2 . = 1. The Hamitonian H2 in terms of the 
Bogoliubov quasiparticles reads 



H 2 = J2 E ^ C W 



(13) 
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FIG. 2: One magnon dispersion. Top: linear spin-wave dispersion 
Z?k of from Eqn. J 14b - Middle and bottom: real and imaginary part 
Re(Im)El i of one magnon dispersion to 0(1/5) from a solution of 
Eqns. (10) - (12) of Refpii on a lattice of 252 x 252 k-points with 
artificial line broadening of r\ — 0.05. 

and the dispersion is given by 

= y/A* -Bl = v/(l - + 2v k ) . (14) 

The magnon dispersion is depicted in Fig. [2] It vanishes 
at the center of the zone, k x — 0, k y = 0, where v^ = 1 and 
at the corners of the BZ, where v^ = —1/2. There are two 
van-Hove singularities, i.e. at E = 3/(2\/2) — 1.061 from 
the maximum energy and at E = 2/3 ~ 0.6667 J from the 
zone-boundary. 

To treat the interaction between magnons we need to ex- 
press the triplic and quartic part of the Hamiltonian, H3 and 
i?4, in terms of the quasiparticles using the transformation 
of Eqn. (fTTT i. For the triplic part we obtain 

H 3 = J2 (ck44- P /(k,P) + 
k,p 

c kC P c -k+ P g (k, p) + 0(c f c f c f + /i.e.)) . (15) 

Terms with three creation(destruction) operators are present in 
principle, but are not expressed explicitly for notational sim- 
plicity. As will be come clear in section [IV] they play no role 



in evaluating the magnon interactions within the Raman re- 
sponse. Moreover 

/(k, p) = iV3(^ p (M k Uk-p + Vk«k-p)(u P + v p ) - 

^(^k + Uk)(%>Dk-p + Uk- p V p ) + 

^k-p(wk-p + Uk-p)(«k% + VkVp)) , (16) 
g(k, p) = i\/3(Pp(Mk-pWk + Uk«k-p)(u P + v p ) + 

^k-p(wk-p + Wk-p)(wkWp + V^Vp) - 
^(Wk + Wk)(Mk- p U p + VU-pVp)) . (17) 

For the quartic part we obtain 

Hi = ~ ~[&S ^ c k c -k c P c - P 

k,P 

+0(c < c 1 c 1 c 1 + c 1 c 1 c 1 c +h.c.) , (18) 

where again terms irrelevant for the Raman scattering are not 
displayed explicitly and 

h(k, p) = 2((t£t£ + vlvl)(vk + 4^ k _ p + vp) - 3(ul+ 
vi)u p v p (2vk + Up) - 3(u p + Wp)w k Wk(^k + 2u p ) 

+4u k v k u p u p (2 + v k + v p + 2f k+p )) (19) 

Eqns. ( fT3T > - ([T9T l allow to construct all vertices relevant to 
the final state two-magnon interactions in the Raman scatter- 
ing. Apart from that Eqns. © - ( [TUt can be used to derive 
the one-magnon selfenergy to 0(1/ S), This has been done in 
Ref. Oil to which we refer the reader for details. For the pur- 
pose of the present work it is sufficient to employ Eqns. (10) 
- (12) from Ref tH to calculate the renormalized magnon dis- 
persion E£ to 0(1/ S). Fig. |2] (middle and bottom panel) 
shows the result of such calculations. It is evident that in 
the real part of the magnon energy, the interactions lead to 
extended and almost fiat regions with a shallow 'roton'-like 
minimum along the BZ faces. Moreover, as the ImE£ al- 
most vanishes at these regions, the life time of a quasi-particle 
with corresponding momenta is very large. On the other hand 
quasi-particles with near-maximum energies are located in the 
momentum regions of rather large damping. 

III. LOUDON-FLEURY VERTEX 

We use the framework of the Loudon-Fleury (LF) model 
for the interaction of light with spin degrees of freedom for 
the calculation of the two-magnon Raman scattering. The LF 
vertex has the form 

R = Eutfin ■ 5){e ont ■ S)SA+s , (20) 

where the polarizations e- m = cos 6x + sm 9y and e out = 
cos 4>x + sin <fiy of the incoming and the outgoing light are 
determined by angles 9 and (j>, defined with respect to the 
x-axis. To derive the final form of the scattering LF vertex, 
we first write spin operators in terms of Holstein-Primakoff 
bosonic a-operators ||3}, and then express the latter in terms 
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FIG. 3: Diagrams for the Raman intensity: a) Bare Raman vertex R 
from Eqn. J2U ; b) Raman susceptibility (both bare, G°, and dressed, 
G, magnon propagators will be considered (see text), c) The integral 
equation for the dressed Raman vertex T in terms of the irreducible 
magnon particle-particle (IPP) vertex 7. d) Leading order 1/S con- 
tributions to the IPP vertex. 



of the boson quasi-particle operators c. We get the following 
expression 



R=J2k M k(ckC-k + c k c 

where Mk is given by 



k c -k> 



M k = (Fi 

3 
4 



F 2 (k,0,0))u k t; k - 



(21) 



(22) 



and we have introduced the following notations: 



F 2 (k,6,cf>) = 25(/ 3 (M)cosfc,+ 
h(0, 0) cos(4f + 3@ky) + h{B, 4>) cos(4f - &ky)) , 
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In principle the Raman vertex contains also c^Ck terms. How- 
ever, at zero momentum and to lowest order in 1/5 these 
terms do not contribute to the Raman response at finite fre- 
quency, and we dropped them. Note that R is explicitly Her- 
mitian. 



IV. RAMAN INTENSITY 

We now calculate the Raman intensity including one- and 
two-magnon renormalizations up to 0{ 1/S). The Raman in- 
tensity 1(H) is obtained via Fermi's golden-rule from 

I{uj n ) = const x Im[ / dr e lUmT {T t {R{t)R))] (24) 



by analytic continuation of the bosonic Matsubara frequencies 
u) m = 1-nmT onto the real axis as iuj m — > O + in, where 
£1 = uji n — uj ou t refers to the inelastic energy transfer by the 
photon, and for the remainder of this paper we assume the 
temperature T = 1/(3 to be zero. The prefactor 'const 7 refers 
to some arbitrary units by which the observed intensities are 
scaled. 

The role of interactions is summarized in Fig. [3] Two 
effects have to be distinguished, namely renormalizations of 
the one-magnon propagators, i.e. G° — ► G, and vertex cor- 
rections to the Raman intensity (final state interactions), i.e 

i?->r. 

All propagators are expressed in terms c^-type of Bogoli- 
ubov particles. To orders higher than 0((l/5) ) the prop- 
agators of these particles are not diagonal i.e. both normal 
G cc (k,r) = -(r r (c k (r)c k )) and anomalous L» cc (k,r) = 
— (?V(ck(T)c-k)) propagators do occur. However, anoma- 
lous propagators are smaller by one factor of 1/5 as com- 
pared to the normal propagators and, therefore, can be ne- 
glected. To first order in 1/5, the normal propagators read 
G(k, iu) n ) = l/(iu> n — ££), i.e. the quasi-particle residue 
remains unity, and E£ is taken from Eqns. (10) - (12) of Ref. 

In order to evaluate Raman intensity (Eqn. d24l i) we have to 
calculate the Raman susceptibility (T t (R(t)R)). In principle, 
the latter can contain terms of type (T T (r _ (r)r - )), with r^ 1 
specified in Eqn. ( f2Tb and Fig. [3] a). However, these terms 
need at least one 0(1/ S) interaction-event to occur or require 
anomalous propagators, i.e., they are smaller by one order of 
1/5 and will be dropped. In the following we consider only 
(T T (r + (r)r~)) + (T T (r~ (r)r + )). By Hermitian conjugation, 
it is sufficient to calculate J(r) = (T T (r~ (r)r + )), which is 
depicted in Fig. [3]b). 

Fig. [3]b) shows the two-particle reducible Raman vertex 
r(k, aj n ,oj m ), which includes a series of magnon-magnon in- 
teraction events. It satisfies the Bethe-Salpeter equation ex- 
pressed in terms of the two-particle irreducible vertex 7, de- 
picted in Fig. [3]c) 



r(k,w„,w m )=r (k) + ^7(k,p 



G( P )G(-p,-w )r(p, ) (25) 

In this work we consider only the leading order contribu- 
tions in 1/5 to 7. They are shown in Fig. [3]d). The quartic 
vertex 74 (k, p) is identical to the two-particle-two-hole con- 
tribution from H4 of Eqns.dT8ll-dT9t and reads 



74(k,p) 



2S 



h(k, P ) 



(26) 
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The two addends forming the irreducible vertex 
73 (k, p, u> n , uj ) are assembled from H3 and one inter- 
mediate propagator, and can be written as 

7 3 (k,p,w„,w ) = -^ E (/(k,p)#(-k,-p)x 



k, P 



G°(k - p, iuj - iu> n )c w c_ k c p c_ p 



(27) 



ff(k, p)/(-k, -p)G°(p - k, iu> n - icj ))c k c_ k c p cL p 



where the functions /(k, p) and g(k, p) obey the symmetry 
relation /[g](-k, -q) = -/[ff](k, q). 

To keep 73 to leading order in 1/5 we retain only the zeroth 
order propagators G° for each intermediate line. In principle, 
H3 allows for an additional two-particle irreducible graph, 
with the incoming(outgoing) legs placed into the particle- 
particle(hole-hole) channel and one intermediate line at zero 
momentum and frequency. However, we verified that these 
contributions vanish exactly. 

Due to 73, Eqn. ( f25T > is an integral equation with respect 
to both, momentum and frequency. This is the first major dif- 
ference to Raman scattering from collinear HAFs, where only 
74 exists at 0(1/5). To proceed, further approximations have 
to be made. Here we simplify 73 by assuming the dominant 
contribution from the frequency summations to result from 
the mass-shell of the propagators in the intermediate particle- 
particle reducible sections of Fig. [3]:) 



VjJ. 

-itu, 



n 



(28) 



This approximation will be best for sharp magnon lines and 
the transferred frequencies iuj m close to the van-Hove singu- 
larities of 2£'k. 

In this approximation for 73, the two-particle irreducible 
vertex 7 simplifies to 



7(k, p, uj n , w ) w 7(k,p) = 7 3 (k,p) + 74 (k,p) 
1 f 2£ k _ p /(k,p)g(k,p) 



(29) 



25 ^ (E k - E V Y - El_ v 



Mk,p) 



Now we can perform the frequency summation over uj q on the 
right hand side of Eqn. (l25l ) as well as the analytic continua- 
tion iui m — > Q + ir\ = z. With this T in the latter equation 
turns into a function of p and z only, leading to 



^L k , p (z)r p (z)=r-(k) 
p 

7(k, p) 



z-2E, 



(r) 



(30) 
(31) 



which is an integral equation with respect to momentum only. 
In the rest of the paper the superscript V refers to the case 
when renormalized propagators with E p are taken in the two- 
particle reducible part of the Raman intensity, while E p cor- 
responds to the usage of bare propagators. 

Close inspection of the vertex 7(k, p) shows, that it does 
not separate into a finite sum of products of lattice harmonics 
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FIG. 4: Raman intensity neglecting final state interactions, i.e. re- 
placing F by M in eqn. (33). Scattering geometry: (f> = 9 = 0. 
Number of k-points: N x N. Dashed (BB): bubble using bare 
magnon energies Ek from eqn. dl4t . shown in top panel of fig. [2] 
Imaginary broadening 77 = 0.003 choosen such as to retain visible 
but small finite-size oscillations. Solid (RB): bubble using renormal- 
ized magnon energies to 0(1/ S) obtained from eqns. (10) - (12) 
of ref.— and shown in the middle and lower panel of fig. [2] Finite- 
size oscillations are suppressed by Im[E^\. The absolute scale of 
J(f2) is set to unity, but the relative scale of BB and RB is kept. 



of the triangular lattice. Therefore, Eqns. (f30b-(f3Tb cannot be 
solved algebraically in terms of a finite number of scattering 
channels, but require a numerical solution. On finite lattices 
this can be done by treating Eqn. d30b as a linear equation 
for r p (z) at fixed z. This marks another significant difference 
between Raman scattering from collinear and non-collinear 
antiferromagnets . 

Finally, the expression for the Raman intensity from Eqn. 
can be written as 



/(O) = const x (J(f2) 

k ft + irj — 



2E, 



(r)J 



(32) 
(33) 



We now discuss the Raman intensity for several levels of 
approximations. First, we neglect final state interactions and 
set r k (z) — > Mk. Fig. [4] shows the Raman intensity as a 
function of the transferred photon frequencies O for this case. 
This figure contrasts the Raman bubble with bare propaga- 
tors against that with renormalized ones. Such results can be 
obtained on fairly large lattices, since they do not involve a 
solution of the integral equation ( [301 . but only a calculation 
of the one magnon self-energyii We keep the shift in off the 
real axis deliberately small in this figure, in order to discrimi- 
nate its effect from that of the actual life-time broadening due 
to the imaginary part of ££. 

First we would like to note that we find the line shape to be 
insensitive to the scattering geometry. This is in a sharp con- 
trast to Raman scattering from the square lattice HAF, where 
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FIG. 5: Effect of final state interactions on Raman intensity. Scatter- 
ing geometry: (f> — 8 — 0. Number of k-points: N x N. Imaginary 
shift off real axis is 77 = 0.03. Dashed line (RB): replacing T by Al 
in Eqn. d33t and using bubble with renormalized magnon energies 
El to 0(1/5) obtained from Eqns. (10) - (12) of Retf. Solid line 
(RBV): using dressed vertex T obtained from Eqn. J30b in Eqn. J331 > 
and renormalized magnon energies E^. The absolute scale of I(fi) 
is set to unity, but the relative scale of RB and RBV is kept. 



Raman amplitudes in A\ g , B\ g and B2 g symmetries are very 
different. 

In case of the bare Raman bubble, one can see two well- 
defined peaks, one at energy SI = 3/v2 and one at SI = 4/3 
- both in units of 3 J/2. These energies correspond to 2 times 
that of the maximum and of the BZ-boundary saddle-point 
of the classical spin- wave spectrum E^. Clearly, the domi- 
nant spectral weight does not stem from the k-points at the 
upper cut-off of the linear spin-wave energy but from the BZ- 
boundary. This does not reflect the bare two-magnon density 
of states but is an effect of the Raman matrix element Mu, 
which samples the BZ regions preferentially. 

Switching on 1/S corrections, two modifications of the in- 
tensity occur. First, both maxima are shifted downwards by 
a factor of ~ 0.7 due to the corresponding renormalizations 
of the one magnon energies. Second, as the BZ-boundary 
saddle-point of E^ has turned into a flat region, occupying 
substantial parts of the BZ for ££, the intensity of the lower 
energy peak is strongly enhanced due to the very large den- 
sity of one-magnon states. Equally important, the imaginary 
part of E 7 k is finite in the BZ-region which corresponds to the 
maximal one-magnon energies. This smears the peak at the 
upper frequency cut-off in /(SI) almost completely - as can 
seen from the solid line in Fig. |4] In contrast to that, Im[E^} 
almost vanishes in the flat regions on the BZ-boundary due 
to phase-space constraints 18 , leading to a further relative en- 
hancement of the intensity from there. 

Next we turn to final-state interactions. In Fig|5] we com- 
pare /(SI) from the Raman bubble obtained with propagators 
renormalized to 0(1/ S) and only bare Raman vertices to the 
intensity obtained by including also the dressed Raman vertex 



Fk{z) from Eqn. (l30b . The numerical solution of the latter 
equation requires some comments. Since the kernel Lic !P (z) 
is not sparse and has rank iV 2 x N 2 , already moderate lat- 
tice sizes lead to rather large dimensions and storage require- 
ments for the linear solver. We have chosen N = 69, lead- 
ing to a 4761 x4761 system which we have solved 200 times 
to account for 200 frequencies in the interval SI € [0,2.5]. 
The kernel has points and lines of singular behavior in (k, p)- 
space, which stem from the singularities of the Bogoliubov 
factors u[u](k) in /[g][/i](k, p) and from the energy denomi- 
nators in Eqn. d29l i. In principle, such regions are of measure 
zero with respect to the complete (k, p) -space, yet we have no 
clear understanding of their impact on a solution of Eqn. d30b 
as TV — > 00. In our case, i.e. a finite lattice, we have chosen to 
regularize these points and lines by cutting off eventual singu- 
larities in ik, p . The comparatively small system sizes require 
a rather larger artificial broadering ir\ in order to achieve ac- 
ceptably smooth line shapes. This can be seen by contrasting 
the dashed curve in Fig. [5] and solid curve in FigJU which 
correspond to identical quantities, however for different finite 
systems, 252 x 252 vs 69 x 69. 

The main message put forward by Fig. [5] is that the final 
state interactions lead to a flattening of the peak from the two- 
magnon density of states, transforming it to a rather broad Ra- 
man continuum. This can be understood, at least qualitatively, 
from the RPA-like functional form of the Bethe-Salpeter equa- 
tion. Discarding momentum dependencies and iterating the 
two-magnon bubble times the irreducible vertex 7, leads to 
a renormalization of the intensity by a factor, roughly of the 
form ~ 1/(1 — 7-p(Sl)), where p(Sl) refers to the two-magnon 
bubble. Directly at the peak position of the Raman bubble this 
renormalization factor may get small, thereby suppressing the 
over-all intensity. While exactly the same mechanism is at 
work also for the square lattice HAF, its impact on the spec- 
trum is complete different. In the latter case the peak-intensity 
without final state interactions is at the upper cut-off of the 
Raman intensity. Suppression of this peak-intensity simply 
shifts the maximum intensity to lower frequencies within the 
Raman spectrum. This shift is then interpreted in terms of a 
two-magnon binding energy. Such reasoning cannot be pur- 
sued in the present case. 



V. CONCLUSION AND DISCUSSION 

To summarize, we have investigated magnetic Raman scat- 
tering from the two-dimensional triangular Heisenberg an- 
tiferromagnet considering various levels of approximation 
within a controlled l/S'-expansion. Our study has revealed 
several key differences as compared to the well-known mag- 
netic Raman scattering from the planar square lattice spin- 1/2 
antiferromagnet. 

First, we found that the intensity profile is insensitive to 
the in-plane scattering geometry of the incoming and outgoing 
light at 0(1/ S). This has to be contrasted against the clear 
difference between Ai s and Bi s .2g symmetry for the square- 
lattice case. 

Second, on the level of linear spin-wave theory, we showed 



7 



that the Raman intensity has two van-Hove singularities. The 
less intensive peak is located at the upper edge of the two- 
magnon density of states and stems from twice the maximum 
of the one-magnon energy. This is similar to the square lat- 
tice case. However, the dominant peak is located approxi- 
mately in the center of the two-magnon density of states. This 
peak stems from the Loudon-Fleury Raman-vertex strongly 
selecting the Brillouin zone boundary regions where the one- 
magnon dispersion on the triangular lattice has an additional 
weak van-Hove singularity. This is absent on the square lat- 
tice. 

Next, we calculated the Raman intensity with the one- 
magnon spectrum, renormalized to 0(1/ S), however neglect- 
ing final-state interactions within the Raman process. In this 
case we have obtained a sharp and almost ^-functional Raman 
peak at energy ~ 3 J/2. At this energy the real part of the 
renormalized one-magnon dispersion shows a large plateau- 
region at the Brillouin zone boundary with a roton-like shal- 
low minimum. Moreover, due to phase-space constraints the 
one-magnon life-time is large in this region. Therefore, the 
two-magnon density of states in this region is strongly en- 
hanced, as compared to the linear spin-wave result. In con- 
trast to that, the intensity at the upper edge of the spectrum is 
suppressed further, since the 0(1/ S) corrections lead to the 
significant one-magnon damping. Finally the overall width of 
the spectrum is reduced by a factor of approximately 0.7. 

In a last step, we considered the impact of the final-state 
interactions to 0(1/ S). Due to the non-collinear ordering 
on the triangular lattice, and in sharp contrast to the square- 
lattice case, we find, that even to lowest order the two-magnon 
scattering is neither instantaneous in time, nor separable in 



momentum space. Our solution of the corresponding Bethe- 
Salpeter equation reveals a broad continuum-like Raman pro- 
file which results from a smearing of the intensity of the 
two-roton peak by virtue of repeated two-magnon scattering. 
While, at this order in 1/S the over-all form of the Raman 
profile is reminiscent of that on the square-lattice, one has to 
keep in mind, that in the latter case the position of the maxi- 
mum in the center of the Raman continuum has to interpreted 
rather differently, namely in terms of a two-magnon binding 
effect. 

In conclusion we hope that our theoretical investigation 
will stimulate further experimental analysis of triangular, 
and more generally frustrated magnetic systems by Raman 
scattering. Several novel materials with triangular structure 
have been investigated thoroughly over the last few years, 
among them the cobaltites, Na^CoO^, and the spatially 
anisotropic triangular antiferromagents CS2CUCI4— and n- 
(BEDT-TTF)2Cu2(CN)3^i. To our knowledge however, mag- 
netic Raman scattering on such systems remains a rather open 
issue. 
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